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ABSTRACT 

We present high-resolution cosmological hydrodynamic simulations of three galaxy clusters employing a 
two-temperature model for the intracluster medium. We show that electron temperatures in cluster outskirts are 
significantly lower than the mean gas temperature, because Coulomb collisions are insufficient to keep electrons 
and ions in thermal equilibrium. This deviation is larger in more massive and less relaxed systems, ranging 
from 5% in relaxed clusters to 30% for clusters undergoing major mergers. The presence of non-equilibrium 
electrons leads to significant suppression of the SZE signal at large cluster-centric radius. The suppression of 
the electron pressure also leads to an underestimate of the hydrostatic mass. Merger-driven, internal shocks 
may also generate significant populations of non-equilibrium electrons in the cluster core, leading to a 5% bias 
on the integrated SZ mass proxy during cluster mergers. 
Subject headings: galaxies: clusters: general - intergalactic medium 



1. INTRODUCTION 

In current theories of cosmological structure formation, the 
cold dark matter and baryon fluid undergoes gravitational col- 
lapse, leading to virialized filaments and roughly spherical 
"halos" on a variety of scales. Collisionless shocks play a 
critical role in this process by converting the kinetic energy of 
in-falling baryonic material to thermal energy. This heating 
leads to large reservoirs of T = 10 6 - 10 8 K gas that fills and 
surrounds massive structures, directly tracing the formation 
and evolution of cosmic structure. 

Observational probes of this gas are generally sensitive only 
to the electron component of the plasma. Typically, theoret- 
ical studies assume that the electrons are in thermal equilib- 
rium with the surrounding ions. However, this is not a good 
approximation in the low-density outskirts of galaxy clusters 
due to the extended timescale for electrons to reach equilib- 
rium via Co ulomb collisions ( |Fox & Lo eb 1997; [Ettori &| 
Fabian 1998 ). In the pre-shock intergalactic medium, the bulk 



of the kinetic energy is carried by the heavier ions. Since the 
electrons and ions are only weakly coupled, this leads to an 
electron temperature in the post-shock plasma that is signifi- 
cantly lower than that of ions. The post-shock electrons and 
ions gradually relax to a single equilibrium temperature. 

The low emissivity of gas at densities of 10" 4 - 10" 5 cm" 3 
has heretofore limited X-ray spectroscopic temperature mea- 
surements to regions where electrons are expected to have 
reached equilibrium (Markevitch et al.|1998| [Vikhlinin et al.| 
|2005[|Pratt et al.|2007| ). Advances in the sensitivity ol X-ray 
instruments are pushing the limit of ICM temperature mea- 
surements to a sig nificant fraction of the virial radius (e.g., 
Geor ge et al.||2009| ). Such measurements remain expensive 



with the current generation of instruments, however. 

In recent years the thermal Sunyaev Zel'dovich effect (SZE) 
has been suggested as a possible probe of gas in thi s phase 
(e.g., |Afshordi et al. 12007 HHallman et al^2007||2009) , due to 
its weaker density dependence. Observations using the next 
generation of SZ experiments are ongoing, and are expected 



to provide accurate measurements of the electron pressure be- 
yond the virial radius. Large samples of massive clusters with 
profile measurements beyond the virial radius should there- 
fore become available in the near future. 

In this work we consider the impact of the two-temperature 
structure of the intracluster medium (ICM) on the thermal 
SZE. We employ high resolution simulations of clusters 
which improve on previous studies ([ Takizawa 1999; Chieze 
[et~aLl[T998l |Courty & Alimi| |2004[ |Yoshida et al'| |2005 ; 
Yoshikawa & Sasaki 2006} by combining the efficient shock- 
capturing properties of Eulerian adaptive mesh refinement 
with high mass and spatial resolution such that both the large- 
scale accretion shocks and internal, merger-driven shocks are 
well captured. We show that the lower electron temperature in 
the cluster outskirts leads to a significant underestimate of the 
gas pressure when derived through the SZE. We also demon- 
strate that the two-temperature structure of the ICM is sensi- 
tive to recent accretion history through internal merger-driven 
shocks and discuss the relative importance of these shocks to 
the SZE. 

2. TWO TEMPERATURE COSMOLOGICAL CLUSTER 
SIMULATIONS 

The simulations were pe rformed using the ART code 
(Kra vtsov et aLl[T997l |2002) , which was modified to model 
a two temp erature electron-ion plasma. Following Yoshida 
|et al.| ( [2005] ), we divide the plasma into two components, elec- 
trons and ions, which are assumed to be individually in local 
thermodynamic equilibrium (LTE) with separate Maxwellian 
velocity distributions defined by temperatures T e and T t , re- 
spectively. The separation into only two species, electrons and 
ions, each in separate equilibrium is reasonable, since the self- 
equilibration timescales for each specie, tu and t ee , are consid- 
era bly shorter than the electron-ion equilibration timescale, 
t ei ( |Fox & Loebll997|[SpTtzer|l962t . 1 When an electron-ion 
plasma passes through a shock, moj 
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TABLE 1 
List of Simulated Clusters 



Name 1 


^200 

(/r 1 Mpc) 


M 2 oo 
(lO 14 //- 1 Mpc) 


T 2 
1 m 

(keV) 


Dynamical state 3 
(Rel/Unrel) 


CL101 .. 


1.77 


12.83 


6.40 


Unrelaxed 


CL104 .. 


1.42 


6.69 


5.44 


Relaxed 


CL6 


0.95 


1.98 


2.14 


Relaxed 



1 Cluster labels correspond to those used in Nagai et al. ( 2007a 1. 

2 T m is the gas-mass weighted temperature measured within r2oo- 

3 Classification of dynamical state is described in Nagai et al. ( 2007b 1. 

goes into heating the heavier ions, causing T t » T e . After 
the shock, electrons and ions slowly equilibrate via Coulomb 
interactions, each converging to the mean gas temperature, 
Tgas = (n e T e + riiTi)/(n e + n{), over a typical electron-ion equi- 
libration timescale, t ei . The evolution of the electron temper- 
ature is given by, 



dT e 
dt 



■ (y-l)T; (V-v), 



(1) 



where the second term accounts for adiabatic compression 
heating and cooling. The timesc ale for equipa rtition between 
two charged species is given by Spitzer ( 1962), 
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where m, T, and Z are the mass, temperature, and charge 
of each specie, respectively, n is the number density, and 
In A « 40 is the Coulomb logarithm. For the fully ionized 
ICM, including contributions from both protons and He ++ , the 
timescale for equilibration is, 



t ei ~ 6.3 x 10 8 yr 



(t;/io 7 k) 



3/2 



(^/10" 5 cm" 3 ) (In A/40)' 



(3) 



Note that this timescale can be comparable to the Hubble time 
in regions with T ~ 10 7 K and overdensities 10 - 100 with 
respect to the cosmic mean. 

While our model assumes negligible electron heating 
within shocks, some non-adiabatic heating due to plasma in- 
stabilities is likely. However, theoretical expectations for the 
amount and sourc e of electron heatin g in these shocks vary 
widely (see, e.g., Byk ov et al.||2008[ for a recent review), 
ranging from a constan tTraction T e ~ 0.27/ of the pre- shock 
kinetic energy (Cargill & Papadopoulos 1988) to a constant 
post-shock electron temperature a nd inverse- square scaling 



with shock velocity T e /Tt ~ v s 2 ([Ghavamian et al. 



2007). 



In either case we expect the electron heating in the high Mach 
number accretion shocks to be small. There are some indi- 
cations of more rapid electron equilibration in sh ocks with 
lower Mach number in the ICM ( Markevitch 2006 ), however, 
the process is not yet well constrained by observations to date. 
We therefore neglect it in our model, noting that our results 
represent a maximum of the possible effect. 

In this work we re-simulat e three galaxy cluste rs selected 
from the sample presented in Nagai et al. (2007a, N07 here- 
after). The initial conditions and simulation parameters are 
identical to N07, except the simulations we present here ne- 
glect the physics of radiative cooling and star formation, 
which should have negligible impact in the cluster outskirts. 
The cluster properties at z = are summarized in Table [T] 
These clusters are chosen to cover a range of mass and the 
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Fig. 1. — The distribution of the mean gas temperature (r gas ) and the ratio of 
electron and mean gas temperatures (T e /T gas ) for three simulated clusters at 
z = 0. Left panel: the projected mass-weighted gas temperature in \h~ l Mpc 
slices centered on each cluster and 12/T 1 Mpc on a side. Right panel: the 
ratio of the electron-temperature T e to the mean ICM temperature r gas . 

dynamical histories. CL101 is a massive, dynamically ac- 
tive cluster, which experiences violent mergers at z ~ 0.1 
and z ~ 0.25. CL104 is a similarly massive cluster, but with 
a more quiescent history. This cluster has not experienced 
a significant merger for the past 6 Gyrs, making it one of 
the most relaxed systems in the N07 sample. The last clus- 
ter in our sample, CL6, is a fairly relaxed, low-mass cluster 
at z = 0, but experiences a nearly equal mass major merger 
at z = 0.6. Each cluster is simulated using a 128 3 uniform 
grid with 8 levels of refinement. CL101 and CL104 are se- 
lected from 120/T 1 Mpc computational boxes and CL6 is se- 
lected from a 80/T 1 Mpc box, achieving peak spatial resolu- 
tion of « 3.6/T 1 kpc and « 2.4/T 1 kpc, respectively. The dark 
matter particle mass in the region surrounding the cluster is 
9 x 10 8 hr x M Q for CL101 and CL104 and 3 x 10 8 h~ l M Q for 
CL6, while the rest of the simulation volume is followed with 
lower mass and spatial resolution. 

3. RESULTS 

Figure [T] illustrates the complex and highly aspherical dis- 
tributions of gas temperature (left panel) and non-equilibrium 
electrons (right panel) in the outskirts of the simulated galaxy 
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Fig. 2. — Profiles of ICM electron temperature relative to the mean ICM gas 
temperature for each cluster at at z = 0. Top panel: 3D electron temperature 
profiles plotted averaged in spherical shells and scaled to r2oo- Bottom Panel: 
2D mass-weighted electron temperature profiles averaged in cylindrical an- 
nuli projected 60/T 1 Mpc through the simulation volume. This quantity is 
equivalent to the bias in the observed SZ flux in each annulus with respect to 
the flux one would observe if the electron temperature, T e , were equal to the 
mean gas temperature, r gas . 

clusters at z = 0. These maps highlight a complicated net- 
work of non-equilibrium electrons associated with accreting 
and shocked material. The recent mergers in CL101 generate 
both internal and external shocks, and these shocks are re- 
sponsible for the significant quantity of non-equilibrium elec- 
trons throughout the cluster. Non-equilibrium electrons are 
found mostly in the cluster outskirts for the relaxed systems 
(CL104 and CL6), and electrons and ions in the inner regions 
are mostly in equilibrium. 

The top panel of Figure [2] shows the spherically averaged 
3D radial profile of T e /T gas as a function of the cluster-centric 
radius for the simulated clusters at z = 0. For the relaxed 
clusters, electrons and ions are in nearly perfect thermal equi- 
librium within r = 0.5r2oo- Beyond this radius, the electron 
temperature becomes increasingly smaller than the mean gas 
temperature, and the magnitude of the bias is inversely pro- 
portional to the ICM temperature: T e is biased low by 5% at 
y - ^200 for CL104 and r = 2r2oo for CL6. The bias reaches 
its maximum at the virial shock, r S h 0C k = [5 - 6]r 2 oo- Beyond 
the shock radius, electrons and ions in the unshocked gas are 
largely in equilibrium, with a slight departure from equilib- 
rium due to the accretion shocks of large-scale filaments. 

The magnitude of the non-equipartition between electrons 
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Fig. 3. — The bias in the hydrostatic mass estimate of galaxy clusters due 
to the non-equipartition of electrons and ions. The bias is defined as the 
difference in the hydrostatic masses computed using the electron temperature 
and the mean gas temperature. 

and ions depends sensitively on the cluster dynamical state. 
The largest deviation from thermal equilibrium is seen in the 
most massive and least relaxed system (CL101). For this clus- 
ter, we find T e /T gas « 0.75 - 0.9 at relatively small radii 
r = [0.5 - l]r2oo and reaching a minimum T e /T m « 0.35 
at r = 3r 2 oo- The non-equilibrium electrons at small radius 
are primarily generated by internal shocks associated with the 
recent merger (see also Figure [T]). 

Production of non-equilibrium electrons by internal shocks 
within the virialized regions of clusters is illustrated in Fig- 
ure [4] which shows the evolution of T e /T gas for CL6 during 
the period immediately following a nearly equal mass merger 
around z = 0.6 (see also Figures 1 and 2 in Nagai & Kravtsov 
2003). The final coalescence generates a strong internal shock 
which propagates out from the cluster center exciting a thin, 
trailing layer of non-equilibrium electrons. Although at fixed 
radius this is a transient phenomenon, lasting for only about 
10 7 - 10 8 yrs due to the short equilibration timescale within the 
high-density ICM, the shock propagation occurs on a dynam- 
ical timescale. This leads to non-equilibrium electrons in the 
cluster center for 0.5-1 Gyr during the merger. These internal 
shocks can therefore introduce a significant temperature bias 
on considerably smaller scales than the accretion shocks. 

The suppression of the electron temperature also leads to a 
bias in the hydrostatic mass estimate, Mhse(< r) oc -dP/dr, 
as shown in Figure [3] We follow the procedure described 



Lau et al.] ( |2009| ) for computing the hydrostatic mass pro- 



m 

files. Briefly^ we reduce spurious fluctuations in the numerical 
derivatives by removing gas in bound substructures and then 
smooth the resulting profile with a Savitzky-Golay filter. At 
small scales, r < 0.7 r2oo, the deviation from the hydrostatic 
mass estimated using r gas is small for all three clusters, and 
remains < 5% for the relaxed clusters out to 1.5 r 2 oo- The 
bias in the hydrostatic mass is quite large for clusters under- 
going recent major mergers (e.g., CL101), where the hydro- 
static mass is underestimated by 20 - 30% in the radial range 
[1-2] r2oo- Note that this effect is in addition to the bias intro- 
duced by the n on-thermal pressure due to random gas motions 
dLau et al.|2009] ). 

These biases in temperature due to non-equilibrium elec- 
trons have important implications for the interpretation of 
Sunyaev-Zel'dovich effect (SZE) observations. This is illus- 
trated in the bottom panel of Figure [2j which shows that the 



4 



RUDD & NAGAI 



z = 0.53 


z = 0.48 


z = 0.43 


z = 0.38 


^rft ***** 









n 1.0 

0.8 
0.6 
0.4 
0.2 



Fig. 4. — The evolution of T e /T gSLS surrounding CL 6 from z = 0.53 to z = 0.38. This cluster experiences a major merger at z ~ 0.6, which is nearly in the plane 
of the figure. This merger causes symmetric merger shocks which propagate outward from the center, eventually reaching and merging with the cluster virial 
shock at z ~ 0.3. 



presence of non-equilibrium electrons leads to significant sup- 
pression of the SZE signal at large cluster-centric radius. We 
define the bias to be the ratio of the Compton-y parameters, 
defined as ysz(T) = J n e Tdl, computed alternately using T e 
and r gas . Again the bias is largest for the unrelaxed system, 
and is as large as 20% in the radial range [0.7 - 2]r 2 oo- The 
bias is more moderate (5 - 10%) for the other two relaxed 
systems in a similar radial range. 

The integrated Compton-y parameter, F S z (or SZ flux), is 
also affected by the non-equilibrium electrons. But, the bi- 
ases are much smaller since the effect is averaged over the 
entire cluster and is dominated by the central regions where 
electrons are in equilibrium. Examining the evolution of the 
bias in 7sz within the spherical region of rsoo for all three sim- 
ulated clusters, we find the bias is no larger than 6%, which 
occurs during periods of major mergers (e.g., CL101 since 
the last major merger at z ~ 0.1 or CL6 at z ~ 0.6 as shown 
in Fig [4]). The bias reaches 9% during the same period when 
measured within the larger radius, r2oo- This process may also 
affect the scatter in the mass-observable relations. At various 
times during the last CL101 merger, however, the cluster lies 
both above and below the mean Y-M relation, while the bias 
due to non-equilibrium electrons always acts to suppress the 
SZ signal. Large samples of simulated clusters are therefore 
needed to adequately address this effect. For comparison, we 
computed the bias in X-ray spec troscopic temperature using 
the spectroscopic-like formula of Mazzott a et al.|(2004| ). The 
resulting bias is negligible, even during mergers, since it is 
strongly weighted to gas at higher densities and correspond- 
ingly snorter t ei . 

4. CONCLUSIONS 

We use simulations of cosmological cluster formation to ex- 
plore the two temperature structure of the ICM and its effects 
on the Sunyaev-Zel'dovich effect. We show that electron tem- 
peratures are lower than the ion (or mean gas) temperature in 
the low-density outskirts of galaxy clusters, where Coulomb 
collisions are insufficient to keep electrons and ions in ther- 
mal equilibrium. This leads to a decrease in the SZE relative 



to predictions which assume electron-ion equilibrium. The 
suppression of electron pressure in turn leads to an underesti- 
mate of the hydrostatic mass. 

Our simulations also show that the magnitude of the non- 
equipartition between electrons and ions depends sensitively 
on the cluster dynamical state. The electron temperature 
(or pressure) is smaller than the mean gas temperature by 
about 5% at r = r2oo for relaxed clusters, but the bias could 
be as large as 30% for unrelaxed systems, in which popu- 
lations of non-equilibrium electrons are created within the 
virial regions of cluster by internal, merger-driven shocks. Al- 
though these are short-lived in the high-density ICM, these 
non-equilibrium electrons can affect global properties, such as 
the integrated Compton-y or the mass weighted temperature, 
by about 5%. X-ray spectroscopic temperature, on the other 
hand, is not significantly affected even during the merger pe- 
riods, since it is strongly weighted to gas at higher densities 
and correspondingly shorter t ei . These biases are the sources 
of systematic uncertainties that will need to be taken into ac- 
count for the interpretation of upcoming large-scale SZE ex- 
periments, including ACT 2 , Planck 3 , and SPT 4 . 

We lastly note that the details of electron equilibratio n will 
affect the abil i ty to detect virial s hocks using the SZ ( |Koc- 
|sis et al.|2 005 ; Mol nar et al.|2009| ) by smoothing the electron 
temperature profile in the region of the shock. Successful de- 
tections of accretion shocks by ALMA 5 can be used to place 
strong constraints on the otherwise uncertain physics operat- 
ing at these shocks. We will explore this possibility in an 
upcoming paper. 
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